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£^ ' Abstract 

a '. 

In this paper we are interested in the barotropic compressible Navier-Stokes sys- 
^vq ■ tern endowed with a non-local capillarity tensor depending on a small parameter e 

such that it heuristically tends to the local Korteweg system. After giving some 
physical motivations related to the theory of non-classical shocks (see [28]) we prove 
0^ global well-posedness (in the whole space M. d with d > 2) for the non-local model 

and we also prove the convergence, as e goes to zero, to the solution of the local 
Korteweg system. 
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1 Introduction 



1.1 Presentation of the models 



This introduction is widely inspired by the works of F.Coquel, D.Diehl, C.Merkle, and 
C.Rohde and we refer to [32], [9] for a deep presentation of the capillary models. 

The mathematical description of liquid- vapour phase interfaces has a long history and 
has been recently renewed in the 80's after the works of Dunn and Serrin (see [14]). The 
first investigations begin with the Young-Laplace theory which claims that the phases 
are separated by a hypersurface and that the jump in the pressure across the hypersur- 
face is proportional to the curvature of the hypersurface. The main difficulty consists in 
describing the location and the movement of the interfaces. 

• I— i 

Another major problem is to understand whether the interface behaves as a disconti- 
nuity in the state space (sharp interface, SI) or whether the phase boundary corresponds 
to a more regular transition (diffuse interface, DI). The DI ansatz is often favored as 
it covers topological changes in the phase distribution (such as the separation of bub- 
bles/drops). Moreover, instead of one system per phase coupled with free-boundary 
problems required by the SI model, in the DI approach only one set of equations has to 
be solved in a single spatial domain (the density takes into account the different phases) 
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which considerably simplifies the mathematical and numerical study. 

Another approach corresponds to determine equilibrium solutions which classically 
consists in the minimization of the free energy functional. Unfortunately this minimiza- 
tion problem has an infinity of solutions, and many of them are physically wrong (some 
details are given later) . In order to overcome this difficulty, Van der Waals in the XlX-th 
century was the first to add a term related to the surface energy to select the physically 
correct solutions, modulo the introduction of a diffuse interface. This theory is widely 
accepted as a thermodynamically consistent model for equilibria. 

The study of the dynamics of a liquid-vapour mixture is complicated: the model is 
based upon the compressible Navier-Stokes equations with a Van der Waals state law for 
ideal fluids. Endowed with such a state law, the compressible Navier-Stokes equation can 
describe the dynamics of multi-phases (liquid and vapour) but the effect of surface tension 
is still not taken into account. In this goal, Korteweg extended (in [26]) the Navier-Stokes 
equations by adding a capillarity tensor modeling the behavior at the interfaces. The 
additional term is directly related to the surface energy added to the free functional in or- 
der to obtain physically relevant minimizers. This local Navier-Stokes-Korteweg system 
belongs to the class of diffuse interface models and, to penalize the density fluctuations, it 
involves the density gradient and in fact it introduces, in the classical model, third-order 
derivatives of the density (and therefore numerical complications). 

Alternatively, another way to penalize the high density variations consists in applying 
a zero order but non-local operator to the density gradient [30], [31]. Let us mention in 
passing that it was the original idea of Van der Waals in [37] and that statistical mechanics 
tend to show that this non-local approach is in fact the correct one (see [32]). This model 
does not introduce additional orders of differentiation but a non-local integral operator 
which creates comparable numerical difficulties. For new developments we also refer to 
the work of M. Heida and J. Malek in [24]. 

In this introduction we will briefly recall both of these approaches and explain their main 
differences. 

1.2 Basic modelization setting 

We want here to explain and give the main lines of how were modeled the local and 
non-local Navier-Stokes-Korteweg systems. Let us consider a single-component fluid 
cointained in an open set Q € R d at constant temperature T*. In the isothermal setting 
the density is denoted by p : f2 — >]0, b[ (b > 0), and the pressure is given by p : p g]0, b[*-> 
p(p) G]0, +oo[. To denote the two phases, we need to define the energy density function 
as: 

W(p) = pf(p,T*), pe]0,b[, and /eC 2 (0,6). (1.1) 
From the mathematical point of view, W has to satisfy the following properties: 

• 3ai,«2 e]0, b[: a.\ < 012, W" > in ]0, ot\ [U]a2, b[, W" < in ]«i,a 2 [- 

• lim^o W(p) = lim p ^ b W(p) = +oo. 
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• W > in ]0,b[. 

We refer to [30] and [9] for a far more complete presentation (and also figures). Since we 
consider the fluid as a single system, the density will play the role of an order parameter 
in the sense that it will determine the phase state. More precisely the high variation of 
density will define the location of the interfaces and we will say that the fluid at the posi- 
tion x is in the vapour/elliptic (or spinodal) /liquid phase if p{x) g]0, ai]/[ai, 02]/ [0:2, b[. 
According standard thermodynamical theory, we define the pressure as follows: 

P(p) = pW'(p) - W(p). (1.2) 

Thus the pressure is a non-convex and non-monotone function ( because P (p) = pW" (p)). 
In passing, let us mention that the free energy / cannot be directly measured while the 
pressure can be determined from experiments. A standard model for the pressure in 
two-phase mixtures is the well-known Van der Waals law given by: 

where R is the specific gas constant and a, b are positive constants. Let us briefly say 
that the parameter a controls the attractive forces between molecules of the fluid and b 
is related to the molecule size. 

1.3 Variational Approach for sharp interfaces 

In the sequel, we want to model our fluid by considering solutions as static equilibria 
at constant temperature T* (at least in an asymptotic configuration in time). In this 
optic, we can associate to this problem of equilibrium for the mixture of liquid-vapor a 
variational problem that consists in minimizing the free energy functional. For this, we 
introduce an admissibility set for the density which takes into account the fact that our 
interfaces have null thickness, more precisely we set: 

A = {p G L?-(£l)/W{p) G L 1 ^), / p(x)dx = m}. 

Jn 

Here we have prescribed the total mass by a constant m > 0. To find a static equilibrium, 
we need to minimize, when p G A , the functional: 

F°[p] = [ W(p(x))dx. (1.4) 
Jn 

The problem is that it is possible to build an infinite number of minimizers: if f3\ G]0, a\[ 
and 02 G]«2;/3[ are the Maxwell states (the two points where W is minimal), up to a 
translation we can have W{j3\) = W(/?2) = and then if < m < /^l^l every 

piecewise constant function p taking its values in {^1, ^2}, so that the total mass is m, is 
a minimizer. This allows drops and bubbles (which are physically correct solutions: they 
only take the values /3i and fa, and away from the boundary, the length of the phase 
interface is minimal) as well as an infinity of minimizers describing an arbitrary large 
number of phase changes (provided that the total mass is m) which clearly are physically 
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wrong. 



In other words, we need selection principles to identify the physically relevant solu- 
tions. 

Remark 1 All the previous minimizers realize the phase transition as a jump between 
the vapour or liquid zones, without taking values in the spinodal zone (at least up to set 
of zero Lebesgue measure). This is why we talk about sharp-interface model. 

1.4 The local diffuse interface approach 

Van der Waals seems to be the first who tried to overcome the lack of uniqueness for the 
sharp interface approach by adding a term of capillarity (see [37]). His idea for selecting 
some relevant physical solution is to penalize the high variations of density (typically 
at the interfaces) by adding some derivative terms to the previous functional. More 
precisely, he considers the admissibility set A local by: 



We now search a function p £ which is a minimizer of the local Van der Waals functional: 



where 7 > is a capillarity coefficient and e > is a scaling parameter. Obviously, func- 
tions which exhibit jumps across a surface are not in H 1 ^). By consequence, a possible 
two-phase minimizer p e cannot change of phase without passing by values in the elliptic 
region, another way to express this phenomena is to say that the constituent propagates 
continuously in the interfaces. This is the reason why the variational problem is called 
a diffuse interface approach. Here the thickness of the interface is expected to be of size 



This diffuse interface approach may appear artificial at first sight but a result, ob- 
tained by Modica in [29], validates this method. He proves that if a sequence of mini- 
mizers p £ converges when e tends to zero, the limit is a physically relevant minimizer of 
the sharp-interface functional F°. 

1.5 The non-local diffuse interface approach 

This alternative approach (also called global diffuse interface approach, introduced by 
Serrin et al in [14] and next by Coquel, et al. in [9], and Rohde in [31] and [32]) also 
consists in a modification of the sharp-interface functional Fq. One of the main interest 
of this approach is that it does not need to introduce high derivatives on the density to 
penalize the high fluctuations of density. 
Let us consider the following admissibility set: 



A iocai = H i(fL)nA°. 




0(e)- 



A9 iobai = A°nL 2 (n), 
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Let us choose a function 4> £ L l (JsL d ) such that: 

(|.| + |.| 2 )<K-) € L 1 ^), / 0(x)dx = l, ^ even, and > 0. 

Such a function is called an interaction potential, and (f) £ (x) = jz4>{^ ) I s called scaled 
interaction potential. Then for 7 > 0, we search for p £ £ A global that minimizes the 
following modified functionnal called non-local Van der Waals functional: 

^ io y[ft] = (w(Pe(x)) + J ^ e (x - y)(p £ (y) - p^fdy^dx. 

Roughly speaking, the non-local term penalizes high fluctuations on an e-scale if the 
interaction potential has most of its mass in a ball centered around zero and whose 
radius is of order e. The set A 9lobal contains functions with jumps but it can be proved 
that minimizers of F^y take values also in the elliptic regions: the model belongs to 
the class of diffuse-interface models. The counterpart of the result of Modica was proved 
by Alberti and Bellettini ([1]), which also validates this approach. 

Remark 2 Let us consider F^ lobal with ft = R and assume that p € C°° (R) is an analytic 
function with the expansion (we refer to [32]): 

1 d k 
p{y) = p{x) + ^ (y - xf—^pix), 



k>l 



If we put this expression into the non-local term in F^ lobal and perform the change of 



variable r = we obtain: 



I I </»(r)(V 7T (er) fc ^ T p(x)) 2 drdx« f f s 2 <P{r){r^-p{x)) 2 drdx. 
Jr Jr. ^[ k\ dx k J R J R dx 

Since J R tfi(r)r 2 dr < +00 holds by definition of <j), the last expression is equal to the local 
penalty terms in Ff ocal (up to a multiplicative constant). In particular the scaling with 
respect to the small parameter e fits. 



1.6 Non classical shocks and the Euler problem 

Another important approach using the capillarity concerns the discontinuous solutions 
of nonlinear hyperbolic systems or conservation laws. Indeed the mathematical model 
of liquid-vapor flows should have special solutions that can be interpreted as dynamical 
phase transitions. In [32], the author shows that there are traveling- wave solutions for 
one-dimensional versions of the local and non local NSK system that connects states in 
different phase. This fact is one of the major arguments to accept the NSK systems as 
promising candidates to model the dynamics of liquid vapour flow in a reliable way. 

We also refer to the works of [3] in the case of the local diffuse approach where the 
authors prove the existence of planar traveling waves representing either diffuse interfaces 
or solitons. The system of equations governing planar traveling waves reduces to a planar 
Hamiltonian system, for which a portrait analysis exhibits heteroclinic /homoclinic orbits 
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(see for more details [3]). 



The first order model that governs the sharp-interface limit is of particular interest 
since it is not a standard conservation law (we will explain in the sequel more precisely 
what it means). In particular it is interesting in this sense to consider the solutions of 
Euler systems as limit solutions of the NSK systems. 

In the sequel, we will interest in the Euler or sharp-interface case, we look for a specific 
volume r in (|, oo) and a velocity v that satisfy in (0, +00) x R the equations: 

d t T - d x v = 0, 

1.5 

d t v-d x (P(r)) = 0, 
with the function P : (^, 00) — > (0, 00) given by: 

P{r) = P{\), TG(loo). 

More precisely if we investigate the well-known Euler system with a Van der Waals 
pressure state in one dimension, we can recall that following the theory of conservation 
law; the two eigenvalues of the system are: 

\i(t,v) = ->J-P'(t), \ 2 ( t ,v) = -^-P'(t). (1.6) 
The corresponding eigenvectors r±, r 2 are: 

= ( fW^) ) ' W2(T - V) = ( -sf^¥\f) ) (L7) 

Furthermore by calculus we obtain: 

P"(t) -P"(t) 
V\ 1 {t,v)-w 1 {t,v) = V ' , V\ 2 {t,v)-w 2 (t,v) = V ; (1.8) 



-P'(r) 2^-P'(r) 

We now recall the definition of a standard conservation law in the sense of Lax (it means 
entropy solutions): 

• The system is strictly entropic if the eigenvalues are distinct and real. 

• The characteristics fields are genuinely nonlinear if we have for all (r, v), 

VAi(r, v) • u>i(t, v ) / and VA2(r, v ) • u>2(t, v ) 7^ 0, 



for more details we refer to [35]. The definition of genuine nonlinearity is some kind 
of extension of the notion of convexity to vector-valued functions (in particular when 
we consider the specific case of the traveling waves). The previous assumptions aim at 
assuring the existence and the uniqueness of the Riemann problem ( see [15] and [35]). 
When P is a Van der Waals pressure, we observe that then the first conservation law 
([35, 15]) is far from being a standard hyperbolic system: 
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• It is not hyperbolic (but elliptic) in (^-, ^) xl, 

• the characteristic fields are not genuinely nonlinear in the hyperbolic part of the 
state space. 

Here the classical theory cannot be applied. It has been shown (see [15]) that there 
exists an infinity of weak entropy solutions of the Riemann problem with initial states in 
different phases. It means that the entropy inequalities are not sufficiently discriminating 
when one characteristic field fails to be genuinely nonlinear. One of the main research 
in this field is to understand what are the physical relevant solutions when 
the flux is non genuinely non linear. An idea for selecting some relevant solutions 
of the Euler system with a Van der Waals pressure is to consider the solutions limit of 
solutions of NSK systems when the capillarity and the viscosity coefficients tends to 0. 
This condition is called the viscosity-capillarity criterium. By a study of traveling waves 
(see [32] or [28] in the context of the Korteweg de Vries system) a remarkable fact is to 
observe that generally the limit solutions violate the Oleinik conditions, more precisely 
the shock are undercompressive (it means non classic in the sense of the Lax theory). 
We refer for more details to the book [28] of P. Lefloch. 

An important research line (see [31], [32], [9]) is to model the capillarity tensor and to 
understand how fast the solutions converges to the Euler system when the capillarity 
and the viscosity coefficients tends to zero. We want point out here that it exists three 
different regimes, more precisely if we assume the viscosity coefficient equal to e with 
e — > 0. Then we have the three different regimes: 

• k << e 2 , the viscosity dominates. 

• k ~ e 2 , intermediary regime. 

• k » e 2 , the capillarity dominates. 

A second important problem from a numerical point of view is to describe the link between 
the local and the non local Korteweg system and more particularly how it plays on the 
choice of the thickness of the interfaces. We propose us to answer to this second question 
by giving a speed of convergence between the local and non local Korteweg system after 
suitable choices of kernel for the local system. 

1.7 Results 

In this work we wish to unify the two models, the local Korteweg system and the non 
local Korteweg system, more precisely we will search to approach the local Korteweg 
solutions by solutions of the non local Korteweg system in an appropriate way (it means 
by choosing suitable kernel on the non local tensor of capillarity). 

Let p and u denote the density and the velocity of a compressible viscous fluid. As 
usual, p is a non-negative function and u is a vector-valued function defined on M. d . In 
the sequel we will denote by A the following diffusion operator 

Au = pAu + (A + ^)Vdiv u, with p > and v = A + 2p > 0. 
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The Navier-Stokes equation for compressible fluids endowed with internal non-local cap- 
illarity introduced in [31] reads: 



(NSR) 



j d t p + div (pu) = 0, 

1 d t {pu) + div (pu <8> u) — Au + V(P(p)) = KpV(4> * p — p), 
supplemented by the following conditions on the kernel 4>: 

(\.\+\.\ 2 )4>(.) ^ L l (WL d ), I (f>(x)dx = l, even, and > 0. 
The mass, momentum equations for the Korteweg system read ( see [14]): 



(NSK) 



d t p + div (pu) = 0, 

d t (pu) + div (pu ®u) -Au + V(P(p)) = upVAp, 



Both of these systems have been studied in the context of existence of strong solutions 
in critical spaces for the scaling of the equations. For example, concerning the strong 
solutions, we refer to [13] for a study of (NSK) system, and to [17] for (NSR). 

If we compute the Fourier transform in the capillarity terms (simplified by the den- 
sity), we obtain (</>(£) — l)p(£) in the non-local model, and — |£| 2 p(£) in the local model. 
We want to see if, when formally <p(£) is "close" to 1 — |£| 2 , we can expect the solutions 
of these models to be close. 

The aim of this paper is to approximate the local Korteweg model (NSK) with a 
non-local model such as system (NSR) introduced on his modern form by C. Rohde in 
[31] and also [9] (see Van der Waaals [37] for the originally works) . For that we will 
choose a specific function (f> £ in the capillarity tensor, and the system we will consider in 
this paper is the following: 

( d t pe + div (p £ u £ ) = 0, 
dt(peUe) + div (p £ u <g> u £ ) - Au e + V(P(p £ )) = p £ ^V((p £ * p £ - p £ ), 



where we set: 

1 , ,x 



-t^(-) with (j)(x) 



£ d ^V ^ ; (2vr) d 

As the Fourier transform of <p is <f>(^) = e - '^' 2 we have: 

and for a fixed £, when e is small, ^ £ ^\ — - is close to — 1£| 2 . 

Remark 3 We mention here that the choice of <f> £ is in accordance with the physical 
relevant capillarity coefficient (see [32]). 

We will consider a density which is close to an equilibrium state ~p and we will introduce 
the change of function p = p(l + q). By simplicity we take p = 1. The previous systems 
become: 



(K) 



dtq + u.Vq + (1 + q)di\ u = 0, 

d t u + u.Vu - Au + P'(l).Vq - nVAq = K(q).Vq - I(q)Au, 
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and 

d t q £ + u £ .Vq £ + (1 + q £ )divu £ = 0, 
d t u £ + u £ .Vu £ - Au £ + P'(l).Vq £ - ^V(0 £ *q £ - q £ ) 

= K(q £ ).Vq £ - I(q £ )Au £ , 

where K and I are real-valued functions denned on R given by: 

Remark 4 In the sequel we will sometimes rewrite the term K(q).Vq as V(G(q)) where 
G is a primitive of K. 

Remark 5 When p ^ 1 the only changes are in the viscosity and capillarity coefficients, 
and in the expression of functions I and K. 

We now want recall differents results concerning the local and non-local Korteweg system 
by emphasizing on the results of global weak solutions. 

Results on the global weak solutions for the local Korteweg system 

The existence of global weak solution for the model of Korteweg with constant capillary 
coefficient for N > 2 is still an open problem. Indeed the main difficulty states in treating 
the capillarity tensor and more specially the quadratic terms in gradient of the density, 
as a matter of fact it seems impossible to pass to the limit in this term as |Vp| 2 is only 
in L 00 ^ 1 ). We can only hope a convergence in the sense of the measure. A way to 
overcome this difficulty would be to obtain regularizing effects on the density. In this 
spirit D. Bresch, B. Desjardins and C-K. Lin in [6] got some global weak solutions for the 
isotherm Korteweg model with some specific viscosity coefficients fi(p) = Cp with C > 
and X(p) = 0. By choosing these specific coefficients they obtain a gain of derivatives 
on the density p where p belongs to L 2 (H 2 ). Recently the second author in [19] has 
shown the existence of global weak solution with small initial data in the energy space 
for specific choices of the capillary coefficients {n{p) = ^ with K > 0) and with general 
viscosity coefficient. Comparing with the results of [6], he gets global weak solutions with 
test function tp G C^ D (]R d ) independant of the density p as it is the case in [6]. In fact he 
is able to obtain a gain of derivative on the density p by using the specific structure of 
the capillary tensor. In particular it is related with the Gross-Pitaevski system via the 
Madelung transformation (see [19] for more details). 

Results on the existence of strong solution for the local Korteweg system 

Let us briefly mention that the existence of strong solutions for iV > 2 is known since the 
works by H. Hattori and D. Li [22, 23]. R. Danchin and B. Desjardins in [13] improve this 
result by working in critical spaces for the scaling of the equations, more precisely the 

d d y d 

initial data (po,poUo) belong to I? 2 2 i x ^21 (^ ne f &c t that 5 2 2 i 1S embedded in L°° plays 
a crucial role to control the vacuum but also for some reason concerning the multiplier 



(Re) { 
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spaces). In [27], M. Kotschote showed the existence of strong solution for the isothermal 
model in bounded domain by using Dorea-Venni Theory and H 00 calculus. In [18], the 
second author generalizes the results of [13] in the case of non isothermal Korteweg system 
with physical coefficients depending on the density and the temperature. He gets strong 

solutions with initial data belonging to the critical spaces i? 2 2 1 x \ x B% \ when the 
physical coefficients depend only on the density. 

Existence of global weak solution and strong solutions for the non-local Ko- 
rteweg system 

The existence of global weak solutions was obtained by the second author in [21] by 
following the Lion's theory. In passing, we would mention that the capillarity allows to 
treat general pressure of the form P(p) = ap 1 with 7 > 1, indeed we get automatically 
some estimate on p in Lf oc what allows us to get renormalized solutions without efforts. 
The first result of well-posedness in finite time are coming from the works of C. Rohde in 
[31]. More recently the second authors in [17] improved these results by obtaining strong 
solutions in initial data invariant for the scaling of the system. 

1.8 Littlewood-Paley theory and Besov spaces 
1.8.1 Littlewood-Paley theory 

In the sequel the Fourier transform of u with respect to the space variable will be denoted 
by F(u) or u. In this section we will state classical definitions and properties concerning 
the homogeneous dyadic decomposition with respect to the Fourier variable. We will 
recall isome classical results and we refer to [4] (Chapter 2) for proofs (and more general 
properties), and to the appendix for specific properties used in this paper. 

To build the Littlewood-Paley decomposition, we need to fix a smooth radial function 
X supported in (for example) the ball B(0, |), equal to 1 in a neighborhood of B(0, |) and 
such that r 1— > x{ r - e r) is nondecreasing over R + . So that if we define <£>(£) = x(£/2)— x(£)> 
then ip is compactly supported in the annulus {£ G M. d , | < |£| < |} and we have that, 

V£GlR d \{0}, 2>(2-'0 = l. (1.9) 

Then we can define the dyadic blocks (A q ) qe z by A q := ip(2~ q D) (that is A q u = 
<p(2~ q £)u(£)) so that, formally, we have 

u = ^A q u (1.10) 

q 

As (1.9) is satisfied for £ / 0, the previous formal equality holds true for tempered 
distributions modulo polynomials. A way to avoid working modulo polynomials is to 
consider the set S' h of tempered distributions u such that 

lim || i S , om||i,oo = 0, 

q^> — 00 
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where S q stands for the low frequency cut-off defined by S q := x(2 q D). If u G S' h , (1.10) 
is true and we can write that S q u = ^ A q u. We can now define the homogeneous 

p<q-l 

Besov spaces used in this article: 

Definition 1 For s G R and 1 < p, r < oo, we set 

.- ^J22 rqs \\A q u\\ r Lp y ifr<oo and \\ u\\g a := sup2 qs \\A q u\\LP- 



\ u \\b s 



We then define the space Bp r as the subset of distributions u G S' h such that \\ 
finite. 

Once more, we refer to [4] (chapter 2) for properties of the inhomogeneous and homoge- 
neous Besov spaces. Among these properties, let us mention: 

• for any p £ [1, oo] we have the following chain of continuous embeddings: 

B° pA ^1/^ B° Pt0O ; 

i 

• if p < oo then B pl is an algebra continuously embedded in the set of continuous 
functions decaying to at infinity; 

• the following real interpolation property is satisfied for all 1 < p,r±,r2,r < oo, 
si / s 2 and 9 G (0, 1): 

i>S2 1 _ f>8s2 + (l-0)s! . 

L-"p,ri ' p,r2l{9,r) p,r ' 

• for any smooth homogeneous of degree m function F on ]R d \{0} the operator F(D) 
maps Bp r in B s v ~ r m . This implies that the gradient operator maps B p r in B^ 1 . 

The following lemma (referred as Bernstein's inequalities) describes the way how deriva- 
tives act on spectrally localized functions. 

Lemma 1 Let < r < R. A constant C exists so that, for any nonnegative integer k, 
any couple (p,q) in [l,oo] 2 with q > p > 1 and any function u of LP, we have for all 
A > 0, 

SuppuC B(0,XR) => \\D k u\\ Lq < C k+1 X k+d( p~^\\u\\ LP ; 
SuppSc {£ G R N /rX < |£| < RX} => Cr* -1 A fc ||ii||£j> < \\D k u\\ LP < C k+1 X k \\u\\ LP . 

The first Bernstein inequality entails the following embedding result: 



Proposition 1 For all s G R, 1 < p\ < P2 < oo and 1 < r\ < r 2 < oo, the space B, 

-<---) 



Pi,ri 



is continuously embedded in the space B P2 ^ P ' ' 1 
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In this paper, we shall mainly work with functions or distributions depending on both 
the time variable t and the space variable x. Most often, these functions will be seen as 
defined on some time interval / and valued in some Banach space X. We shall denote 
by C(I;X) the set of continuous functions on / with values in X. For p £ [l,oo], the 
notation L P (I; X) stands for the set of measurable functions on / with values in X such 
that t i-> \\f{t)\\x belongs to W(I). 

In the case where / = [0, T], the space L p ([0,T]; X) (resp. C([0,T];X) will also be 
denoted by L^{X) (resp. Ct{X)). Finally, if / = E + we shall alternately use the notation 
L p (X). 

One of the main advantages of using the Littlewood-Paley decomposition is that it 
enables us to work with spectrally localized (hence smooth) functions rather than with 
rough objects. As a consequence, we naturally obtain bounds for each dyadic block in 
spaces of type L^{L P ). Going from those type of bounds to estimates in Lj,(B pr ) requires 
to perform a summation in £ r (Z). When doing so however, we do not bound the L^,(B pr ) 
norm for the time integration has been performed before the l r summation. This leads 
to the following notation (after J.-Y. Chemin and N. Lerner in [8]): 

Definition 2 For T > 0, s € E and 1 < r, p < oo, we set 

\\ U \\l p T {B^ r ) := ll 2,7S |l A 9 n llL£,(LP)ILr( Z )- 

One can then define the space lf T (B, ' ) as the set of tempered distributions u over 
(0,T) x R d such that linw__ 00 S„u = in Z/([0, Tl; L°°(E d )) and ||u|| 7 p^ s , < oo. The 
letter T is omitted for functions defined over IR + . We shall also adopt the notation 

Cr{B 9 p , r ) ■= L?(B s p!r )nC([0,T};B s P;r ) and C 6 (B» r ) := L°°(B s p>r ) nC b (R + ; B^). 

The spaces Lj,(Bp T ) may be compared with the spaces Lj,{B p r ) through the Minkowski 
inequality: we have 

W u he T {B lr) < ll«llL5.(B., r ) X r>p and ||«|lz^(^ r) > W u \\ L ^B kr ) if r < P- 

The general principle is that all the properties of continuity for the product and compo- 
sition which are true in Besov spaces (see below) remain true in the above spaces. The 
time exponent just behaves according to Holder's inequality. 

Let us now recall a few nonlinear estimates in Besov spaces. Formally, any product 
of two distributions u and v may be decomposed into 

uv = T u v + T v u + R(u,v) (1-11) 

with: 

T u v := S q -iuA q v, T v u := S q -ivA q u and R(u, v) := A q uA q /v. 

11 1 \q'-q\<\ 

The above operator T is called "paraproduct" whereas R is called "remainder". The 
decomposition (1.11) has been introduced by J.-M. Bony in [5]. 

Let us now briefly recall the following estimates (we refer to [4] section 2.6, [10], [17] 
for more properties of continuity for the paraproduct and remainder operators, sometimes 
adapted to Lf^(Bp r ) spaces). These estimates will be of constant use in the paper. 
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Proposition 2 There exists a constant C > such that for any s £ 1 and any p,r G 
[1, oof, we have for any (it, u) G x B* , 

\\T u v\\ 6 . <C 1+ l s l||«|| L oc||i;||^ . 

p,r P,t 

Moreover, for any s£l, any t < 0, for any p, ri , 7"2 G [1, oo[, we have, for any (it, u) G 

Ot y DS 

"oo,ri A p,r%> 

C 1+ l s+ 'l , 1 <fe/ 11 

||T u u|| r> s +t < \\u\\m \\v\\ns with - = min(l, ! ). 

D P,r — I " D oo,r 1 - D p,r 2 f T\ Tl 

And concerning the homogeneous remainder: 

Proposition 3 There exists a constant C > which satisfies the following inequalities. 



Let si,S2 € I and pi,P2, n, ?"2 G [1, oo]. Let us assume that s\ + S2 < | or s± + S2 < | 
if r = 1, and 



1 de/ 1 1 ^ , , 1 def 1 1 ^ , 

- = 1 < 1 and - = 1 < 1. 

p Pi p 2 r ri r 2 

If si + S2 > then we have, for any (u, v) G Bp 1 iri .B p "^ r2 , 

(j\s 1 +s 2 \+l 

\\R(u,v)\\fys l+ s 2 < ■ \\u\\b°i \\v\\b s 2 ■ 

B p,r g\ + S2 PI - r l " a P2, r 2 

When r = 1 and Si + S2 > 0, we have, for any (u,v) G Bp^ ri .Bp^ r2 , 
\\R(u,v)\\ A s 1+S2 < Cl Sl+S2 l +1 || U ||^ 1 ||«||^ a . 

c p,oo - D Pl,ri" - D P2. r 2 

Using Besov injections together with the previous estimates, we get the following result 
(we refer to [4] or [12]): 

Proposition 4 There exists a constant C > which satisfies the following inequalities. 
For any si, s 2 G R, any 1 < pi,P2,P < oo and any 1 < ri,r 2 ,r < oo such that 

111 ,111 

si + s 2 >0, -< 1 < 1, and -< 1 < 1, 

p p\ P2 r r\ T2 

we have 

Cl Sl+S2 l + 1 def 111 

\\R(u,v)\\tfi,2 < ; IMInn \\v\\a*2 With (71,2 = s 1 + s 2 -d( 1 ), 

" B P,r Si + S 2 "^Pl.V nti P2^2 X p x p 2 p 

provided that a± t 2 < d/p, or &\ 2 < d/p and r = 1. 

In this article we will frequently use the following estimates (given by the previous recalls) 
under the same assumptions there exists constant C > such that: 

II^uHIrs < cildU 00 Hull r S < cm d Hull ™ , 

ii u ii ,o — ii in- ii n±( 2 — ii ii ii \\H 2 y 

D 2,l 

\\T u v\\ os+t < C\\u\\ o< ||u||™ < C||u|| t+ d||u||o S (r < 0), 

II <* ll±( 21 — II II -tfJo oo II M « 2 i — 11 n -0 2 ,l V 7 

' D 2,l 
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\\R(u,v)\\^ 1+ s 2 <C\\u\\^J\ v \\^ <C\\u\\ d \\v\\ 6 '2 (s 1 + s 2 >0), 
\\R(u,v)\\ d < C||«||^i \\v\\gs 2 < C\\u\\ 6 'i \\v\\tf2 (s! + s 2 >0). (1.12) 

R 2 ,oo 2 ,± Z,L 

D 2,l 

Let us now turn to the composition estimates. We refer for example to [4] (Theorem 
2.59, corollary 2.63)): 

Proposition 5 1. Let s > 0, u G £f,i n L °° and F 6 W^ 2 ' 00 ^) such that F(0) = 
0. Then F(u) G an( ^ there exists a function of one variable Cq only depending 
on s, d and F such that 

WF^W^KCoiWuh^WuW^. 

2. Ifu and v G B\ x and if v - u G for s G] - |, f] and G G ^ [ s J +3 '°°(M d ) such 

that G"(0) = 0, then G(u) — G(u) belongs to B| ;1 and there exists a function of two 
variables C only depending on s, d and G such that 

\\G(v) - G(u)\\ 6 . <C(||u|| L oo,||v|| LO o) \\u\\ d + \\v\\ d \\v-u\\g. . 

\ B 2,l B 2,l J 

Remark 6 As 2. is proved using 1., together with the fact that G(v) — G{u) = (v — 
u). Jq G'(u + t(v — u))dr and the previous estimates on the paraproduct and remainder, 
we can prove that, more generally, 

||G( V )-G(u)||^ <C(\\u\\ L ~,\\v\\ L ~) (\G'(0)\ + \\u\\ d +\\v\\ d ) \\v-u\\ 6 . 
1.8.2 Hybrid Besov spaces 

As announced before we will get a different regularity for the density fluctuation in some 
low and high frequencies, separated by a frequency threshold. This will define the hybrid 
Besov spaces. Let us begin with the spaces that are introduced by R. Danchin in [10] or 
[4]: 

Definition 3 For a > 0, r G [0, oo] and s G R we denote 

lei 

In the present paper, we will only use these norms with r G {l,oo}: 

ihi 5 .,oo = Yl 2 ' (s_1) ii^Hl2 + Yl a2ls \\^i u h*, 

;<io g2 (i) ;>io g2 (^) 



and 



i<l°g 2 (i) «>log 2 (i) a 
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Remark 7 As stated in [4] we have the equivalence 



^(iHI^+alHI^J 



< \\U SS.OO < U Ri-1 | Gi I Hi I DS • 

— ll n±( a — H "£> 2 ,1 -^2,1 



And if we denote INI^-inss = IMIbj- 1 + IMIs;^' then 



— min(l, a)||u||^ s -i n g S < ||m||^s,oo < max(l, a)||it|| 



1 

2 

In this article we will mainly use the following hybrid Besov norm: 



Definition 4 Jf l £ is a frequency threshold, and s, t £ R we define the following hybrid 
norms: 

\\u\\^ t d = f £ 2<iAH| i2 + £ ^2"IIMl»- (1-13) 



Remark 8 In the following we wiJJ take Z e = l^og 2 (-^p) — 1] as defined in (2.36). It 
goes to infinity as e goes to zero. 

Remark 9 If e is small enough, we have \\u\\^s,i < ||u||^s,i < ||n||g S +i, s 
1.9 The Korteweg System 

Unlike the Navier-Stokes or (NSR) systems, here the density, as the velocity, is regu- 
larized for every frequency (in fact there is a frequency threshold in the Fourier modes, 
but in both cases the density is parabolically regularized). In the sequel we choose ~p = 1 
and we will denote qo = po — 1 and q = p — 1 the density fluctuations. We will recall 
here some results from [13] about the strong solutions of the Korteweg system (for more 
simplicity, we do not mention here any exterior forcing term): 

Theorem 1 Assume that P'(l) > 0, min(/i, 2/x + A) > 0, that the initial density flue- 

d ^ . d . d 

tuation qo belongs to B 21 Pi i? 2 2 1 , and that the initial velocity uq is in (Bj 2 ! ) d . Then 
there exist constants tjk > and C > depending on k, p, v, P'(l) and d such that if: 

IMI . A- 1 . A + IKH . d_! < TjK 

D 2,l 1 ID 2,1 -°2,1 

then system (NSK) has a unique global solution (p, u) such that the density fluctuation 
and the velocity satisfy: 



1+2x 



q G C(R+, n B^) n L 1 (R + , B*J n B 2 2 7 ), 

nGC(M+,i?|~ 1 ) d nL 1 (R + ,i?|+ 1 ) rf . 



Moreover the norm of (q,u) in this space is estimated by the initial norm C(\\qo\\ <j_ 1 d + 



J 2,l 1 '-"2,1 



l«0|| .$-1, 

■°2,1 



Further in this article R. Dan chin and B. Desjardins provide a Fourier study of the 
linearized system and observe different behaviours wether the quantity v 2 — An is positive, 
negative of zero. In any case they obtain parabolic regularization. 
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1.10 Notations and main results 
Definition 5 The space E§ is the set of functions (q, u) in 

(c b (R+, b^ 1 n B} tl ) n l\r + , B s e +1 ' s n E| +2 > s )) x (c b (R+, b^ 1 ) n l\r + , b s 2 + 

endowed with the norm 

def 

\\(q>v)\\e; = \\u\\ l °°b°-i + IklLoo^-i + lkll L ooB| x 

+ W u \\lib s 2 + 1 + IMI L ib|+M + \\q\\ L ib s +^ s ( 1 - 14 ) 

We first prove global well-posedness for system (NSR £ ) by following similar ideas as in 
[17] and uniform estimates with respect to e: 



Theorem 2 Let e > and assume that min(p, 2p + A) > 0. There exist two positive 

constants ijn an> 

■ --i 

uq G i? 2 2 i and 



2 — 1 . £ 

constants r]R and C only depending on d, k, p, A and such that if qo € B 21 n -B 2 2 1; 



Ikoll . d-i . d + IKII < 

B 2,l nB 2,l B 2,l 

d 

then system (NSR £ ) has a unique global solution (p, u) with (q, u) G E£ such that: 
\\(q,u)\\ d < C(\\q \\ a ! d + ||n || d J. 



Remark 10 Note that in the low frequencies, the parabolic regularization for q is the 
same as for the Korteweg system, that is the low frequencies of q are in B 21 n B 21 . 



The main result in this article is the following: when the initial data are small enough 
(so that we have global solutions for (NSK) and (NSR £ )) the solution of (NSR £ ) goes 
to the solution of (NSK) when e goes to zero. 

, d ^ . d 

Theorem 3 Assume that min(/x, 2p + A) > 0, P'(l) > and that q G B 2 2 x n B^ v 
■ --i 

uq G B 2l . There exists < rj < min(nx, "i?) such that if 

IMI .i-i -4 + IK II .4-1 ^ r ?> 

then systems (NSK) and (NSR £ ) both have global solutions, and with the same nota- 
tions as before, there exists a constant C = C(r\, K,~p, P'(l)) > such that for all a G]0, 1[ 
(ifd = 2) or«G]0,l] (ifd>3), 

ll«e-«||_ . i- a -i + \\9e ~ Q\\~ .4-«,-i + l|9e-9||- . %- a 
L? B h LTBh L t° B 2,i 

+ \\u £ - «||- l6 |_ a+ i + life " 9||_ 1 ^_ a+1 ,^_ a + life - ?|| z ^|_ a+2 ,d_ a < CV\ (1.15) 
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Remark 11 The same results hold for any function <p G S(M. d ) such that V£ G M. d , 
m = g(\Z\ 2 ) with: 

• Function g : R + — > R taices its values in [0,1], with g(0) = 1, 

• Function h : x i-> 1 ~^ a: - > is decreasing with limo h = 1 and lim^ h = 0, 

• Function k : x ^ 1 — g(x) is increasing with limo = and lim^ fc = 1, 

• For all 1 < (3 < 2, there exists Cp > such that for ah x > 0, 

ff(x) - 1 + x 
° " ^ " ^ 

The structure of this article is the following: in the second section we will obtain a 
priori estimates on the linear system with convection terms. The third section is devoted 
to the existence and uniqueness of solutions for the nondocal model, and in the last 
section we will obtain the convergence result of theorem 3. In the appendix, one will find 
the proof of some estimates involving the special hybrid Besov norms introduced in this 
paper. 

2 A priori estimates 

In this section we focus on the following linear system (e > is fixed and for more 
simplicity we write (q,u) instead of (q £ ,u £ )): 

{dtq + v.Vq + div u = F, 
k 
d t u + v.Vu - Au + pVq — -^V((j) £ * q - q) = G. 

With 

Au = fiAu + (A + u)Vdivu. 
This section is devoted to the proof of the following a priori estimates: 

Proposition 6 Let e > 0, s G R, I = [0,T[ or [0,+oo[ andv G L l (I, B^ 1 )nL 2 (I, B^ ± ). 
Assume that (q,u) is a solution of System (LR £ ) defined on I. There exists a constant 
C > depending on d, s, [i, v, p, k, cq and Cq such that for all t G /, 

IMIZfs^ 1 + IMIz-b*- 1 + IMIz^B^ + II'"IIz ( i bj+ 1 + IMIzib £ s+i ' s + Ikllzisf 1 - 2 ' 8 

C7/ *(||Vw(t)|| a +\Hr)\\ 2 d )dr 

<Ce B ^ (||«o||^-i + Ikob-i + WqoWb^ 



+ \\n^^\n nHl + \\G\\ n ^). (2.i6) 
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The proof of the proposition will be close to those in [10], [4], [11] or [17]: lead by the 
behaviour of the linearized system and using symmetrizers we get estimates on dyadic 
blocks and obtain the expected estimates. The difference here is that we need to carefully 
estimate penalized terms and we obtain, for the density, a frequency threshold depending 
on e. 

We choose here not to separate the compressible and incompressible parts, and we 
will localize in frenquency: using the Littlewood-Paley decomposition, for all Z £ Z we 
define qi = Aiq and ui = A/ii. We obtain the following system: 

{d t qi + v.Vqi + div m = F t + R h 
k 
d t ui + v.Vui - Am +pVqi - ^ V(0 e *qi- qi) = Gi + R h 

where Ri = [v.V,A[]q and R[ = [v.V,A[]u. 

If a > is fixed and small enough (for example we can choose a = min( / u, A + 2/x)/4), 
we will introduce in the following: 

_ k 

h t = \\ u l\\h +p\\Qi\\l2 + a H V ^Hl2 + 2(^|V^) L2 ) + -^(qilqi - 4> £ * 

2.1 Preliminary result 

In order to prove the estimate we will show as a first step the following lemma: 

Lemma 2 Under the previous notations, there exist some m > and C > such that 
for all IeZ: 



\jhf + m(2 2l \\ Ul \\ 2 L2 + ||V«||i 2 + ^(Vg,|Vg/ - e * V«)l») 

< C(||V«|Uoo + \\vf L oo)h 2 + C((l + 2 i )(||F i || i2 + WRih*) ( 2 - 17 ) 

+ II^/IIl 2 + \\ R 'i\\L 2 )hi- 



Proof: If we compute the innerproduct in L 2 of the first equation from system (LLR £ ) 
by qi, and the innerproduct of the second equation by u\ we obtain that for all / G Z: 

^H^lli 2 + ( div ^l«)L2 = -(v-Vqi\qi) L * + (Fi + Ri\qi)v, (2-18) 

and 

1 d k _ _ 

'\ u l\\ 2 L* ~ {Aui\ui) L 2 +p(X7qi\ui) L 2 ^(4> £ * Vqi - Vqi\ui) L 2 



2dV 



-{v.V Ul \ Ul ) L 2 + (Gi + RWn)^. (2.19) 



Using an integration by parts allows us to write that (div ui\qi) L 2 = — (Vqi\ui) L 2 and 
then to combine the previous estimates in order to eliminate these terms. We need to get 
rid of terms like this because they cannot be absorbed by the left-hand side and after the 
use of a Gronwall type estimate they introduce multiplicative constant terms e ct (which 
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are problematic as we look for global time estimates). It is less easy to get rid of the 
k _ _ 

term — {4>e * V<# — Vqi\ui) L 2 and we will explain later how to do it. 

Integrations by parts also provide that there exists a constant C > such that: 

\(v.Vqi\qi) L 2\ < C\\Vv\\L°°\\qi\\ 2 L 2, 
\(v.Vui\m) L 2\ < C\\Vv\\l°°\\ui\\ 2 L 2, 
as well as the fact that 

-{Aui\ui) L i = n\\Vui\\ 2 L 2 + (A + /x)||divu;||| 2 

which leads to 

— (Aui\ui) L 2 > v\\X7ui\\ 2 , where v = A + 1\x and v_ = min(/i, v). 

(if A + fi > it is immediate, else we use that ||divitj||^ 2 < II Vu/|| 2 .) Combining these 
estimates we obtain that: 

Id k 

2^dNli 2 + PWllWh) + "\\Vui\\ 2 L 2 - ^(4>e*Vqi-Vqi\ui) L 2 

< C\\Vv\\ L ~(\\uif L 2 +P||«||i2) +P(\\Fl\\ L 2 + ||^|| i2 )|| % || L 2 

+ (||^|| L2 + ||^|| L2 )||^|| L2 . (2.20) 

As we want to obtain regularization on the density, we need at least a term such as ||f//||^ 2 
to appear in the left-hand side. If we compute the inner product of the equation on u\ 
by Vqi we get: 

k 

(d t Ui\Vqi) L 2 + (v.Vui\Vqi) L 2 - {Aui\Vqi) L 2 +p\\Vqif L 2 + ^{Vqi - <t> e * V qi \Vqi) L 2 

= (G l +R[\Vq l ) L 2. (2.21) 

In order to get a full derivative, we need to estimate the term (ui\dtVqi) so we will 
write the equation on Vqi and its innerproduct on uf. 

d t Vqi + V(v.Vgj) + Vdivn; = VF, + VR t . (2.22) 

Before adding (2.21) to 

{d t V qi \ui) L 2 + {V{v.V qi )\ Ul ) L 2 + {Vdiv Ul \ Ul ) L 2 = {VFi + VRi\ Ul ) L 2, (2.23) 

we remark that a simple computation shows there exists a constant C > such that: 

Kw.VujIV®)^ + (V(«.V«)l«i)L a l < C\\Vv\\ L o \\Vq l \\ L 2\\u l \\ L 2. 

Moreover, integrations by part also provide the fact that 

(A Ul \V qi ) L 2 = v( Ul \VA qi ) L 2. 

Gathering these estimates, we obtain that: 

d _ k 

-^(ui\Vqi) L 2 - v(ui\VAqi) L 2 - ||divuj||| 2 +p||V^||| 2 + ^( v 9« ~4>e* Vgj|Vg/) L 2 

< c||v«|| L -||V%|| L2 ||^|| L2 + (||G,|| £a + KMHVall^ + 2'(||F,|| La + WR^WmWv. 

(2.24) 
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k 

Remark 12 Note that —(Vqi — (f> £ * Vg;|Vg;) L 2 is nonnegative. 

In the previous equation, the term (ui\V Aqi) L 2 is annoying because it introduces much 
more derivatives that we will be able to handle in the following, so we need to neutralize 
it. A simple way to do that is to compute the innerproduct of the equation of Vqi by 
Vqi. Using that: 

(Vdiv«,|V«)L2 = NVA%) L2 , 

and that 

\(V(v.V qi )\V qi ) L 2\ < C||V«|| L =o||V«||| a , 
we get the following estimate: 

^j^mWh + HVAm)^ < c||v v ||^||v%|li 2 +2 i (||F/||L2 + ||^|| i2 )||v«|| L 2, (2.25) 

so that, adding v (2.25) to (2.24) will neutralize the previous annoying term. 

Remark 13 Due to — ||divu;||^ 2 appearing in (2.24) we have to be careful and look at 
(2.20)+a (v (2. 25) +(2. 24)) with a > small enough (a < v) so that — a||div«/||| 2 can 
be absorbed by u\\Vui\\ 2 L2 . 

k 

As announced we have now to manage to get rid of — (cf) £ * Vqi — Vqi\ui) L 2 which is 

introduced by estimate (2.20). A rough estimate will eventually provide some constant 

term e~ 2 which will ruin any asymptotic approach as e goes to zero. The easiest way to 

k 

take care of it is to consider the equation satisfied by the nonnegative term -^(qAqi — 
4>e * Qi)l 2 - Computing the innerproduct of the equation on q\ by q\ — (j) e * qi gives: 

^(<#k - <Pe * ~ (ui\Vqi - 4> e * Vqi) L 2 

= -(v.Vqi\qi -(j) £ * qi) L 2 + (F, + Ri\qi - <p £ * qi) L 2. (2.26) 

We need to be careful when estimating the right-hand side as this estimate will be 
multiplied by k/e 2 and added to the others. Using Plancherel (up to a multiplicative 
constant) we introduce the following notation: 

When estimating the right-hand side, we are forced to let the block qi — <p £ * qi in one 
part, as our only possibility is to absorb some terms using remark 12. As the following 
function: 

/ : R* + ^ M 

1 - e~* ( 2 - 27 ) 
xi— > , 

X 

is decreasing from 1 to 0, and as up to loosing some derivatives we want to get uniform 
estimates in e, we can write that for all £ G M. d : 

l_ e - £ W<^| 2> 
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So, using Plancherel, 

\(v.V qi \ qi -4>e*qi)^\ = C\ [ ZVq].y/l - e~^\ 2 Vl - e-e 2 \t\ 2 .W£jdt;\ 

<C [ \^ l \^\.V / l-e-^\\\WM. (2.28) 

Thanks to the fact that we are localized in frequency, and using the Cauchy-Schwarz 
estimate, we obtain that 



\(v.Vqi\qi -4> £ * <7z)l 2 I < Ce2 l \\v.Vqi\\ L 2\\y/l - 4> £ * qi\\\ 2 



< Ce2 , H L oo||V«|| £2 v'(<n|« - * Qi)iP- (2-29) 
Using the same arguments allows us to write: 

m + Rilqi - <Pe*qi) L 'A =C\ I FfTRi.Vl - e^Vl - ^FW.£(f)d£| 

J~R d 

< CetfdlFtWv + \\Ri\\i?W(qi\qi-<l>e*Qi)v>- ( 2 -30) 

Plugging (2.29) and (2.30) into (2.26), we obtain that: 



^{qi\qi-4>e*qi)L2-{uiNqi-4>e*^qi)L2 <Ce IM|L°°||V®|| L 2y 

+ C^dlF,^ + HHI^jySSzJ*^. (2.31) 

Now, when we compute (2.20)+a (i/ (2.25)+(2.24))+ (2.31), we will be able to neutral- 
ize or absorb any annoying term. For that we will introduce as previously announced: 

k 

h l = INIII2 +p\\qi\\h + a { u ¥Jqi\\h + 2 (^|v^)l 2 ) + ^{qi\qi -<t>e*qi)L 2 - 

Using the Cauchy-Schwarz estimate, we can write that 

\(ui\Vqi)\l 2 < II^IIl 2 -I|V^|| L 2 < ^||V^||| 2 + -||«i|||2, 

so 

, 2ck . || .. 9 || n2 olv 1 1 9 k . 1 . 9 

(1 - — )\\ui\\ L 2 +p\\qi\\ L 2 + — II V« 111,2 + -^{qi\qi - 4>e * qi)L 2 < h i 

< (1 + ^)|M|| 2 +p||«||i2 + ^||V«||| a + ^(qi\qi - 4> £ * qi)v- (2-32) 

The last term is nonnegative and lesser than ||V<fo||£,2. If we want hf to be equivalent 
to ||«;||^2 + lk/|| £ 2 + ||Vg/||| 2 we need that a < u/2. Remember that we already have a 
condition on a: for — a||divn;||^ 2 to be absorbed by the left-hand side, we need a < v. 
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So we need a < min^u, |) and for example we will fix here once and for all a = |. After 
computing (2.20)+q (i/ (2.25)+(2.24))+£(2.31) we get: 



\jh 2 i a)||V«,||i 2 + ap\\V qi \\ 2 L2 + ^(Vqi - 4> £ * V qi \V gi ) L 2 

< c\\vv\\ L ~ti} + c ((i + 2 z )(||f / || l2 + n^|| L2 ) + hqii^ + ||i?;|| L2 ) h 



+ Ck2\\ mL , + m \ L ,)^^SMl 



, ^ I, / (Vg;|V« - 4> £ * Vg/) L 2 



Thanks to the definition of hi, we have: 

k^mL* + \\Ri\\l^)\[^^^ < Vk2 l (\\ Fl \\ L2 + 
and using the classical estimate ab < (a 2 + b 2 )/2, we obtain for the last term: 



Cfc||?;||i,oo||Vgj||£,2 



(Vqi\Vqi - (f> £ *Vqi) L 2 
e 2 



ak kC 2 
< 2^(V9i|V«-^ e *Vgi)L» + -2^-||«llL<»||Vgi|||2. (2.34) 



So, if we denote by m = min(z^ — a,ap, ^) > we have proved there exists a constant 
C > such that for all / G Z, 



2~dt 



hf+m (\\Vui\\ 2 L2 + ||V«||£ a + ^(Vqi - & * Vgj|V«)L») 

< C(||V«||l=c + IbUloc)^ 2 + C ((1 + 2')(||Fi|| L2 + \\Ri\\ L 2) + WdWv + ||i?;|| L2 ) hi, 

(2.35) 



which ends the proof of lemma 2. 



2.2 Proof of the proposition 

According to the previous estimate, we observe that each term of hi except Vqi is parabol- 
ically regularized. The key point is that thanks to -^{qi\qi — <p £ * 9/)l 2 (that appears in 
hi) we will be able to get regularization on Vqi for some low frequencies. This leads to 
the hybrid Besov spaces defined in the introduction. 

-The low frequency case: More precisely, the frequency threshold is simply given 
by the study of the decreasing function / introduced in (2.27). As f(z) goes to 1 when 
z goes to zero, there exists some 7 > such that for all < z < 7, f(z) > \. If we 
denote by < cq < Cq the radii of the annulus C used in the dyadic decomposition, let the 
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threshold l £ be the greatest integer such that e 2 2 2 ' e Co < 7 (that is l £ = log 2 (^P^) — 1]). 
Then for all frequency I < l £ and all £ G 2 l C, we have: 

(e|£|) 2 < e 2 C 2 2 21 * < 7 and 1 - > (2.36) 

and then for all I < l E , 

r r 2 r 2 2 21 

(V« - <fc * V«|V^)l 2 = C / (1 - e~ £ l«l )|V«(0| 2 d£ > C— §— ||V«|li a - 

So plugging this into the estimate given by lemma 2 gives: 

+ m ^||v^||| 2 + nvall^ + ^te - 4> £ * qi\qih* + c^\\v qi \\ 2 L2 

< C(\\Vv\\ L ~ + \\v\\ 2 LOO )hf + C ((1 + 2 / )(||F i || i2 + ||i2,|U*) + ||G,|| La + ||^|| L2 ) h h 

(2.37) 

and we have obtained that for all I <l £ : 
S"' + 2 2 

< C(||V»|U» + |H||«)ft, + C((1 + 2')(||fi|U* + + ||G,|U. + KHy) , (2.38) 
that is for all / < l £ and all t: 

hi(t) + m2 21 I hi{r)dT 
Jo 

< h t (Q)+C J (WVvhco + WvWl^ht+C ((1 + tiXmh* + \\Rih*) + IIGiH^ + ||i?||| i2 ) dr. 

(2.39) 

Remark 14 Notice that as l £ ~ l°g2(2c^)' an< ^ £ oes to inanity when e goes to zero. 
The parabolic regularization eventually reaches every frequency. 

-The high frequency case: Remember that l £ is characterized by: 

fe2'«Co < V7. 
\e2^+ 1 C > ^7 ' 

So, for all Z > Z £ + 1 and all £ G 2 Z C, we have 

e 2 |£| 2 > e 2 2 2 'c 2 > e 2 2 2 ^+ 1 )C ,2 (^) 2 > 7 (^) 2 - (2.40) 
Let us go back to estimate 2. We can rewrite it into: 
Id 2 21 1 

--/i 2 + m2 2! (IMH2 + WmWh) + 2^(«l« " & * + ^2 ( v «l V« - * V«)l» 

< C(||Vt;|| L =o + \\v\\ 2 LOO )h 2 + C ((1 + 2')(||F, + ||^|| L2 ) + ||G,|| La + ||i?i|| L2 ) /i,. 

(2.41) 



23 



Using Plancherel, (2.40) together with the fact that function x *— > 1 — e x is increasing, 
we get that: 

(V qi \V qi - cf> £ * V Ql ) L2 = C f (1 - e- £2 l«l 2 )|V^ (0| 2 de > (1 - e- 7( ^ )2 )||V^|li 2 . 

As 2 2/ > ( an d greater than 1 if e is small enough) if we denote by m' = m. min(^j, 1— 

e c o ) > 0, our estimate becomes: for all I > l e + 1, 

^ + 7^ 

< C{\\Vv\\ L oo + Ml^h + C ((1 + 2 / )(||F i || L2 + ||^|| L2 ) + ||C?,|| La + ||^|| L2 ) , (2.42) 
that is for all I > l £ + 1 and all i: 

£ Jo 

h{Q)+C j (||V«||loo + ||«||£oo)/i,+C ((1 + 2 l )(\\Fi\\ L 2 + ||^|| i2 ) + HGjIl^a + ||i?||| L2 ) dr. 

(2.43) 



< 



We need to find a way to get parabolic regularization for U[. Classically (we refer for 
example to [10], [4] or [17]) we will use the new information given by the previous estimate 
in the energy estimate for uf. 

~\\Mb+i^\Mh 

k _ 

< C 1 1 1 1 || || H-( || C/ || ^2 "I- II ^/ II i. 2 ) ||-Lt/ || ||j^2 • ll^ 2 "'"^^ II — ^e:*^^ ||^2 . 1 1 ^ 1 1 2 , 

(2.44) 

and then, simplifying and integrating in time: 
r-t 



\\ Ul {t)\\ L 2 + E 2 21 f \\ Ul {T)\\ L 2dr < |h(0)|| L2 
J o 

+ C jf (\\Vv\\ L ~\\ Ul \\ L 2 + (||G Z || L2 + KH^a + (p+^)\\V qi \\ L2 ^ dr. (2.45) 

Thanks to (2.32), estimate (2.43) implies 

% f \\Vqi{T)\\ L 2dT<C'% fh^dr, 
£ Jo £ Jo 



so that we can estimate the last term of the left-hand side of the previous estimate and 
then there exists a constant such that for all I > l £ + 1 and all t: 

h{t) + v2 21 ( \\ui{T)\\ L 2dT+^ I hi(r)dT 
Jo £ Jo 

< hi(Q)+C J (||V«||l=o + ||«||£oo)/i,+C ((1 + 2')(||F / || L2 + ||^|b) + \\Gi\\ L 2 + ||^|| L2 ) dr. 

(2.46) 
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Remember that: 

hi ~ IKIU 2 + Iklli 2 + l|Vg/|| L 2. 

Let us collect what we have obtained for the low and high frequencies: if we use the 
Gronwall lemma to (2.39) and (2.46) we obtain that for all t, for all I < l £ : 

+ \\<n(t)\\i? + \\Vqi(t)\\i? + ™2 2i f (|h(+>|| £ 2 + \\ qi (r)\\ L 2 + \\V qi{r)\\ L ,) dr 

Jo 

< eC/o'dlv^WlUco+ll^Wllicojdr^^ojii^ + | k/(0) || i2 + || Vffi (0)|| La 

+ C ((1 + 2<)(||^|| £2 + ||i^|| £2 ) + ||C?,|| £a + K|| L2 )dr), (2.47) 

and for I > l £ + 1: 

\\Mm^ + hl{t)\\L^m{t)\\L^+f G (v? 1 \\ui{t)\\v + ^(\\qi(T)\\L* + ||V«(r)|| L2 )) dr 

< eC/o'dlV^rJIUco+ll^rJIlicJdr^^oju^ + || % ( )|| L2 + ||V%(0)|| £ 2 

+ C^* ((1 + 2')(||F i || L2 + ||i^|| £2 ) + ||G,|| £a + ||^|| £2 )dr), (2.48) 

If we use the definition of the hybrid norm introduced in the first section of this paper, 
after multiplying these estimates by 2 S_1 and a summation over / £ Z, we obtain that 
there exists a constant C > such that: 

IMIZf BJ^ 1 + IMIZfBj;- 1 + IMIZfB^ + ll n llz t 1 B=+ 1 + lkllziB| +1 ' s + \\l\\ZjB s e +2 ' s 
< c e c J?(ll Vw M^- + H T )lli-)*(||«,(0)||^-i + 1191(0)11^-1 + IIVwWH^ 

+ 11^11^ + 11^11^+11011^-1+ Aj:(2^- 1 )+^)||i? / |lL 2 +2^ 1) ||^|| £2 )dr). 

Jo zez 

(2.49) 

Recall that i?/ = [-u.V, A/]g and R' t = [v.V, A/]u. We refer for example to [4] (lemma 
2.96 and estimate 10.10) for the following well-known results: 



4+i 



Proposition 7 There exists a constant C > such that for all v £ i? 2 2 1 and g 6 Z>f 1; 
(/+/ = Az]<? £/iere exists a nonnegative sequence c £ / 1 (Z), mi/i q = 1, swc/t £/za£ 



/or a/Z Z 6 



|n|| L 2 < Cc/2 (s ||Vt;|| d llffll^. . 



Then there exists three nonnegative sequences c, c', and c" G t" (Z) (with N c ; =1) 
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such that for all I G Z: 

'||i2i|| L2 <Cqi2-'-||V«|| 4 , 

"2,1 

62 d 2,1 



|itf|| L 2 < C&V^" 1 ) II Vull d ||«||6.-1. 

I illi - i II 11^ II IIB 2jl 



This allows to write: 



<C [ \\Vv(t)\\ .a (\\<i\\ ir , + |M| z} +||»||/r ' ulr ( "- r ' Ui 
Jo 



r2 '-#2,1 " M± >2,1 

J3 1 



Plugging this into estimate (2.49), and using once more the Gronwall lemma finally gives: 

IMIZfBI^ 1 + II^IIl^BJ- 1 + II^IIZoo^* i + + lkllziB|+M + \\q\\ligs+2,s 

Cf*(\\Vv{T)\\ L ^+\Hr)\\loo+\\^v{r)\\ a )dr 

<Ce B ^ (\\u Q \\gs-i 

+ Ikob-i + ||Vgob ri i + + W F hlB h + II G IIlib-0- ^ 2 - 51 ) 

- d 

which, thanks to the injection ^ ^°°! en ds the proof of the proposition. ■ 

3 Proof of theorem 2 

The proof of the global well-posedness of system (NSR £ ) is classical and follows the lines 
of the proof for the compressible Navier-Stokes system (we refer to [13], [10], [4] (section 
10.2.3) or [17]). The only difference here is that we use, when proving the existence, the 
estimate given by Proposition 6, which ensures uniform estimates with respect to e. 

The uniqueness is also classical: when d > 3 we also use estimation (2.16), but when 
d = 2 (as for the Navier-Stokes system) we reach an endpoint for the estimates of the 
remainders in the paradecomposition, and the simplest is to use the same methods as for 
example in [4] (section 10.2.4 and the use of logarithmic estimates). The negative powers 
of e won't cause any trouble as the uniqueness is proved when e > is fixed. 

3.1 Existence 

3.1.1 Step 1: Friedrichs approximation 

In order to use the classical Friedrichs approximation, we define the frequency truncation 
operator J n by: for all n G N and for all g G L 2 ' 
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and we define the following approximated system (we omit the dependency in e for more 
simplicity) : 



dtq-n + Jn {J n U n -^ J n q n ) + J«div U n = F, 



ni 

- y' ( \ T - 7 T HP. 



d t u n + J n (J n u n .X7J n u n ) - AJ n u n + P (p).VJ n q n - -^^{^e * J n q n - J n q n ) = G n , 
where 

[F n = -J n ( J n q n .div J n u n ) 

[G n = J n (K(J n q n ).VJ n q n ) - J n (I(J n q n )AJ n u n ) 

It is easy to check that it is an ordinary differential equation in L n x (L n ) d , where 
L n = {u € L 2 (M. d ), J n u = u}. Then for every n £ N, there exists a unique maximal 
solution in the space C 1 ([0, T* [, L n ) and this system can be rewriten into: 

{d t q n + Jn (u n .Vq n ) + divu n = F n , 
i K P 
d t u n + J n (u n .Vu n ) - Au n + P'(p).Vq n - V(</> £ *q n - q n ) = G n , 

where 

(F n = -J n (q n .divu n ) ^ 
[G n = J n (K(q n ).Vq n ) - J n (I(q n )Au n ) 

3.1.2 Step 2: Uniform estimates 

Using estimate (2.16) with s = ^ from Proposition 6 (as the proof is based on energy 
methods, the operator J n does not have any effect) we obtain that for all t £ [0, T* [: 

C/o'dlVttnMH d +ll«n(r)|| 2 d )dr 

h n (t)<Ce B ^ *h (|K(0)|| + ||g„(0)|| + \\q n (0)|| .* 



I! ■' ' " " " ' " R 2 * " " 1. 7 

JD 2,l D 2,l - D 2,l 

+ 11^11-, .d.i + H^nlL +||Gn||- 1 .d_ A (3.53) 
LiB,?, LiB?, Lis?, / 



where we denote 



. /.\ 11 11 11 11 11 11 

n n {t) = \\u n \\_ . + \\q n \\~ . 4_i+ 9n ~ .4 

L-B^ L t ~B 2 ^ LfB^ 

Let rj > be small (it will be fixed later) and assume that: 

h(0) = f \\u(0)\\. 4_i + l|9(0)|| ,+11,7(0)11 , <r?. 

B 21 B 2,l S 21 

Let us now define 

T n = sup{t G [0,T*[, /i n (t) < 2C7i(0)}. 



27 



As for all n G N, /i„(0) < /i(0) < rj, we have T n > (C > 1) and we will now prove by 
contradiction that T n = T*. Assume that T n < T*, then for all t < T n , h n (t) < 2Ch(Q) < 
2Ci] and: 

/ ||Vu„|| d dr < h n (t) < 2Cr] 

JO 

and by interpolation 

/ \\u n \\ 2 d dr < \\u n \\ . djunll . d +1 < h n (t) 2 < 4C 2 i] 2 . 

JO B*! L°° B li V-Bli 

We obtain that for all t < T, 



Kit) < Ce 2C ^ 1+2C ^ h n (0) + ||F B || d , + \\F n \\ ri d + ||G„|| ^ • 

\ ^ B 2,l ^ S 2,l ^ "2,1 / 

We need now to estimate the right-hand side (see 3.52 for the expressions). Thanks to 
estimates (1.12) and Proposition 5 we obtain that: 

u n\\ . d+i ||(7n|| . d-i an d ||-fn||.d <C||^n||d+i 

"2,1 "2,1 "2,1 "2,1 "2,1 "2,1 

ll^(9n)Vg n || < C\\K{q n )\\ d \\Vq n \\ d, < CJqJ 2 , , 

2,1 B 2,l B 2,l ^2,1 

\\I{q n )Au n \\ d_ , <C\\I(q n)\\ . d \\U n \\ . d +1 n\\ . d \\U n \\ . d +1 . 

^2,1 ^2,1 ^2,1 ^2,1 ^2,1 

So we can write that for all t <T n , 

K(t) < Ce 2c2 ^ 1+2C ^ (h n (0)+ rKH.^dknll.^ + llgnll.- )dr) 

V Jo B li B h B li J 

+ / \\q n \\ 2 d dr. (3.55) 



/ lkn|| 2 d 
JO B^ 

Using proposition 10 (see appendix) we obtain that 

hnf. d < C(\\q n \\ d_! + lkn|| d )\\q„\\ d +h d 

B 2,l 2,1 B 2,l B e 

which gives 

K(t) < Ce 2C ^ 1+2C ^ (h(0) + J\\\u n \\^d +1 + \\q n \\^ +1 d)h n (T)dT 
Then, thanks to the Gronwall lemma and to the fact that 

I (K|| -4+1 + ||9 n || d+i dT ^ M*) ^ 2Ch ^ ^ 2C? ?' 

JO B 2 ,i B ? 
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we obtain that 



/()(!MI <f + 1 + ll<?n|| d +1 d)dr 

h n (t)<Ce 2c2 ^ l+2C ^h(0)e »Zi B? <Cf( V )h(0), 

where 

f(rj) d = e 2c2r ii l + 2Cr i) e 4C 2 V e 2c2 v(WCv) 

This function goes to 1 as rj goes to zero so if 77 > is fixed and chosen so that f(n) < | 
then for all t < T n , h n (t) < 3C7i(0)/2 < 2Ch(Q) which contradicts the fact that T n is 
maximal. Then T n = T* and for all t < T*, we have h n (t) < 2Ch(0) < 2Cr] and in 
particular 

|u„|| . d_ 1 < 00 

L T n B 2,l 

\Qn\\ ^ . d_i < OO 
L T n B 2,l 

And, as J n is a frequency cut-off operator on C(0, ^,n) we have the following bounds: 

f IWIls? l 2 < 00 

/ J-n 

1 lkn||L2? L 2 < 00 

These bounds blow-up when n goes to infinity but all that is important is that for a fixed 
n, it implies (using Cauchy-Lipschitz and the fact that T* is maximal) that T* = 00. So 
for all n the solution (q n ,u n ) is global and we have the uniform bound (in n and e): 

Vt > 0, h n (t) < 2Ch(0) < 2Cr] (3.56) 

Once it is obtained the rest of the method is very classical. We will use the uniform 
estimates to bound the time derivatives of the approximated solutions. Here will appear 
negative powers of e, but as it is fixed, and as only n goes to infinity, the computations 
are strictly the same as in [4] or [10]. Then from this we easily prove that (q n ,u n ) is 

uniformly bounded in C 2 ^^ 2 ' 1 •* x C4( R +> B 2,i ) ( Q f course, uniformly in n but not in s, 
but e is fixed for now) . The Ascoli theorem together with the classical diagonal extraction 
and weak convergence arguments give, for each e > 0, the existence of a global solution 
(q £ ,u e ) whose norm is bounded by 2Cr/: 

Theorem 4 There exists 77 > such that under the same assumptions, if 
IH0)|| . 4_! + ||p(0) - p\\ . i-i + IIp(o) -p\\.i < v 

R 2 R 2 R 2 

D 2,l D 2,l D 2,l 

System (NSR £ ) has a global solution (p, u) such that ifq = p — ~p then (q, u) is a solution 
in E s £ of (Re) and 



||( 9 , «)|| „ 4 < 2C(||g(0)|| j., ^ + ||«(0)|| ,_J < 2Cr?. 

2,1 
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3.2 Uniqueness 



Once again, as e is fixed, the computations are the same as for the compressible Navier- 
Stokes system. As in [4] or [17], we successively use the transport estimate for the density 
fluctuation ([4] theorem 3.14) and the transport-diffusion estimate for the velocity ([4] 
Proposition 10.3). If d > 3, it is simpler to use the a priori estimates provided in the 
present paper. We need to separate the cases d > 3 and d = 2 (the case d = 2 is more 
difficult because of endpoints for the remainder estimates in the Littlewood-Paley theory) 
and we obtain the following result: 

Theorem 5 Let d > 2 and (%, Ui) (i G {1, 2} ) two solutions of (R £ ) with the same initial 

d 

data on the same interval [0, T* [ and both belonging to Ei (T*). Then there exists a > 
such that if, for i G {1,2} 

then (<?i,ui) = (q 2 ,u 2 ) on [0,T*[. 



4 Proof of Theorem 3 

In this section we will show that the solution of (R £ ) goes to the solution of (K) and give 
estimates of the speed of convergence as e goes to zero. For that we will once more use 
Proposition 6. 

As said in the introduction, if the initial data satisfy 

||<7o||.d i . d + 1 1 i*o 1 1 . 4_! < mm(rj K ,m), 
ft 2 nR 2 R 2 

D 2,l 1 l - D 2,l D 2,l 

then systems (K) and (R £ ) both have global solutions (q £ ,u £ ) and (q, u), and with the 
same notations as before, there exists a constant C = C(r], k) > such that for all £ G R 
we have, 

IKIL^ ■ |-i + Ike II ■ |-i + ll&ll-^ ■ i + ll u e|| j-l+i 

L t° B 2,l L t° B 2,l L T B 2,1 L t B 2,l 

+ he\\ L ^i +1 d + he\\ L]§ l + 2d < 2Crj, (4.57) 

and 

1 1 1 1 . d_! + ||<Z||_ . d_i + |k||_ . d +\\u\\ . d +1 

L T B 2,\ L t° B 2,l L t° B 2,l L t B 2,l 



+ lkll Ll ^ + i + Ikll lfi 4 +2 < 2C^. (4.58) 



J t LJ 2.\ ±J t 1J 2,l 



Up to an additional forcing term, let us write system (K) with a capillary term as in 
system (R £ )' 

(Re) 

( d t q £ + u £ .Vq £ + (1 + q £ )divu £ = 0, 

| d t u £ + ti £ .Vji £ - ^ln e + P'(l).Vq £ - ^V(4> £ * q £ - q £ ) = K{q £ ).Vq £ - I(q £ )Au £ , 
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and 

dtq + u.Vq + (1 + q)div u = 



(K) 



d t u + u.Vu -Au + P'(l).Vq - ^V(^ e *q-q) = K(q).Vq - I{q)Au + R £ , 



where R £ d = — ^V(</> £ * q — q — e 2 Aq) and we recall that i<C and I are the following 
real- valued functions defined on R + (whose value is when (7 = 0): 

™ -(*<»- 3^) - '<«>-j*r 

In order to estimate the new term R £ we use the following lemma (which can be proved 
by a simple function study): 

Lemma 3 For all 1 < /3 < 2 there exists a constant Cp > such that for all x > we 
have: 

so that we immediately obtain 

Corollary 1 For all s G R and (3>0,ifq£ B*+ 1+W then R £ G and tiiere exists a 
constant Cp > such that we have 

2,1 £> 2,1 



In the following we will use this estimate for (3 and s such that {5 > 1 and s+2/3 G [|, f+1]- 
Let us now write the system satisfied by the difference (Sq, Su) = (q £ — q,u £ — u): 



d t Sq + u £ .VSq + div Su = 5F, 

d t Su + u £ .X75u - ASu + P'(l).VSq - ^V(4> £ * Sq - Sq) = SG - R £ , 



(4.59) 



where 



with < 









'8G X 


= —SuSJu 


'5F l = 


—Su.Vq 




SG 2 


= (K(q £ ) - K{q)) .Vq £ 


< SF 2 = 


— Sq.div u £ 


and < 


SG 3 


= K(q).VSq 


JF 3 = 


—q.divSu 




5G4 


= (I{q £ )-I{q))Au £ 








JG 5 


= -I(q)ASu. 



We will estimate (Sq, Su) in less regular spaces. Let us show that for a < min(l,<i — 1) 
(a will be precised later): 



' 2 



Sq G LT^lx n LfBl x n L\Bi 
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Like previously, from estimates (4.58) and (4.57) we obtain (using Proposition 10 with 
s = | for q £ , and interpolation for q) that q £ and q are uniformly bounded (with respect 

to t and e) in L/^B^i and that for all p G [l,oo], ii e and it are uniformly bounded in 

1. Classically we use dtq £ = —u £ .Vq £ — divu £ — q £ divu £ , and estimate the right-hand 
side: 

• Taking p = we get div u £ G L\~ a B^ \ ' 

^- • — CK 

• As u £ G L t a -B 2 2 i from (1.12) we get that 
|<? £ divu e || .d_ a _! < C||g e || .d ||divu e || .^_„_ 1 

^2,1 ^2,1 ^2,1 

|« e Vg e || .$-„-i <C||g e || d ||u e || d_ a . 

™2,1 -°2,1 -°2,1 

-z- ■ ^-a-1 

So that Ojg £ is bounded (uniformly with respect to e and t) in L\ a B% 1 and 
then 

q £ G C t ^B|-«-i and ^ _ 9e(0 ) = & - g e LfBll*- 1 
uniformly in e but not in t. As the same is true for q, we get for all t > 

bq^C t Bll a - X nL?Bi; x a -\ 

■ i ■ i~ a ■ <L 

By interpolation, as 5q G L^B^, we also have <5g G C tJ B 2 2 1 n Lf^ 2 x for all t > 

and <fy G L\B* ' 2 D Lj.B e 2 ' 2 (at this stage we cannot get uniform 
bounds in t nor e, and this is not a problem for fixed t and e). 

Remark 15 As we had to estimate \\u £ \\ ,d_ a , this requires | — a G [| — 1, | + 1] 

S 2,l 

that is a € [—1,1]. 

2. Similarly we use 

d t u £ = -u £ .Vu £ + Au £ - P'{p).Vq £ + ^V(0 e * g e - g e ) + K(q £ ).Vq £ - I{q £ )Au £ , 

~2^/>f-«-i 



Taking p = gf- we get „4u e G L^ a B^ 



Remark 16 As we estimate \\Au £ \\ d_ a _ 1 , this requires | — q + 1 G [| — 1, | + 1] 
that is a G [0,2]. 

• As Vq £ G LTBl\ X n LfB^' 2 we have Vg e G L^Bf ~ a_1 . 
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Thanks to (1.12) we prove that 



|u e .Vu e || 2 j} _ a _ 1 < C\\u £ \\ d ||Vu e || 2 (J _ a _ 1 



L t " S 2,l 



||J(g e ).Au e ||_ 2 d_ a _ 1 < CJ&H a \\Au e \\_ 2 . d_ a _ 1 



L t B 2,l 



L t " B 2,l 



\\K(q e )Vq e \\.i_ a _ 1 <C\\K(q l 



r2 

°2,1 



Ml . d- a ||Vg e || d_! 
B.?, B?, 



2,1 



2,1 



We obtain 



l|tf(<fe)V<fe|| .$_„_!< C||K(g £ )|| .$_J|V<fe|| d x . 
Remark 17 As we used Proposition 5 we need | — q > 0. 

~ ■ d„ Q _^ . . 2 _ d_ j 

So dtu £ is bounded in the space LfB^i + £f 2 ~" B 2 2 i (uniformly in t but not 
in e) and then 



-a-l 



u £ e C t 2 S| 1 and u £ - u e (0) = ii e - 6 LfB 2 2 l 



a-l 



Using the same argument used for Sq, we obtain that Su € L^B^i ° 1 r)LjB 2 2 1 



-Q + l 



We can now use Proposition 6 with s = | — a (if a < min(l,d — 1), the following 
computations are true for the cases d = 2 and d > 3), we obtain: 



L t B 2,l 



L T B 2,1 



O/o'dlVueWH d + IMt)|| 2 d )dT 



+ ||^|| £ ^|_ a+1 ,d_ Q + ||(5g|| z; ^d_ a+2j d_ Q < Ce 



2,1 



D t a 2,l ^t a 2,l 

Thanks to (1.12) we estimate the external force terms: 



\5F! 

\SF 2 
\5F 2 
\5F 3 



•$-«-! = ||*«.Vg|| . rf-a-l < C\\q\\ d +1 , \\SU\\ .$-„-!■ 
S 2,l S 2,l S 2,l "2,1 



R2~ 
^2,1 



< c 



■ D 2,l 



R"2" 
D 2,l 



d .11(5^11 d 

R 2~ 11 "oH- 



= H^-diVUeH .d_ a _! < C||^|| d-a-lll^ll,.^!, 
B 2,l °2,1 °2,1 



'b^ 

°2,1 



* a <C\\5q\\ a a \\u e \\ d +1 , 

"1 1 -Do 1 1 



2,1 



2,1 



S 2,l B 2,l ^2,1 ""2,1 
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B 2.1 



• 4 _ a <C\\q\\ - 

"2,1 "2,1 "2,1 

Recall that from Corollary 1, 

B 2,l B 2,l 

And concerning the other forcing terms in the equation on 5u: 

• \\SGiW d_ a _, = \\Su.Vu\\ d_ a _i < C*||<5u|| 

"2,1 -"2,1 "2,1 "2,1 

• ll<SG 4 || = || (/(&) -/(g)M« e || .rf_«_i < C||/(<fe)-/(g)|| . g_J|u e || . 4 

o 2 R 2 r> 2 

"2,1 "2,1 °2,1 

and using Proposition 5, we can estimate the first term: 

\\I(q e )-I(q)\\ <C(||<fe|U~Jg|| LO o) ( |/'(0)| + ||<fe|| | + ||?|| ) ||<fe-?|| . 

B 2,l \ B 2,l "2,1 J B 2,l 

so there exists C v (bounded with respect to rf) so that 
||<5G 4 || a a _ t < C^H d a || 

"2,1 "2,1 "2,1 

• ||«? 5 1| = ||/(g).A<HU-«-i < C\\I(q)\\ a \\6u\\ 

R2 r>2 o 2 r>2 

"2,1 D 2,l D 2,l -"2,1 

and similarly, as 7(0) = 0, || -/"(<?) || . d < Co ( 1 1 ^ 1 1 z.°° ) 1 1 9 1 1 . d so that 

B 2A B 2.1 



\\SG 5 \\ d_ a _ 1 < C v \\q\\ d \\ou\\ d. 

B 2,l B 2,l B 2,l 

The last two terms require more attention: as we want to get global in time estimates, we 
cannot afford to loose any power of t so we need to use any available L^-type estimates 
of q and q £ . 

• If a < d — 1 then we can use Proposition 11 with s = | — a and t = | — 1 (and 
here s + t > 0): 

\\5G 2 1| d_ a _, = || -K(q)).Vq e \\ d_ a _, 

B 2,l B 2,l 

< C (\\K(q e ) - ^(?)ll^|- a -i + II^Gfe) " ^(9)11^-) ll V ^ll^^-i- ( 461 ) 

Moreover as a g]0, d — 1[, we have | — a, | — a — 1 g] — § , f ], and we can use 
Proposition 5 with K'(0) / 0, 



||^(? e )-if(?)||^_ a _ 1 + ||/v(r/J-A'(r/)||_. i4 „ < C„{ \\5c;\\. A , , + INI...4 
and we obtain: 



D 2,l "2,1 \ "2,1 "2,1 



ll^n^.^ < c v ( ||«y g || A 4_ a _ 1 + \\Sq\\^d_ a ] ik e |L, +1 ,,. 



"2,1 \ "2,1 "2,1 / "e 
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• Similarly, Proposition 11 with s = | and t = ^ — a — 1 (s + 1 > 0) implies: 
||*G 3 || , = ll^-V&zH a_ a _, 

B 2,l B 2,l 

< C + II^WII^J llV^II^-a.j.a-x. (4-62) 

One more time, as i^T(O) = 0, Proposition 5 allows to write: 



< C v + , ||<5 9 |l ,_ a+1 ,. 



3 2,1 V S 2"",l °2,1 



Remark 18 As we used Proposition 5 with s = | — 1, the previous estimate is valid 
only when d > 3. 

In the case d = 2, as q G Lf^B^i H L\B^i , thanks to interpolation, we obtain that 
q€L*Bl v 

~ ■ -— a— 1 

Similarly, but instead of using interpolation we use Proposition 10, as <5g € Lf (£ 2 1 n 
n L\Bl~ a+l ^~ a , we obtain 5g G L 2 t Bj~ a . So 

||<5G 3 || .$_„_!< , ||V^|| 

S 2,l °2,1 "2,1 

As||K( g )||.d <C(||g|| L =o)|| g || .d <C„||g||.$ < C||g|| . 4 , we obtain 

B 2,l ^2,1 B 2,l ^2,1 



•°2,1 -B 2 2 ! Kf^ \ #2,1 #2,1 / B £ 2 2 

< ^Ikll .|-llNI .jf-^1,^ + ^Ikll (ll^ll + ll^ll^.-,) • 

Remark 19 As we used Proposition 5 we need the sum of indices s and t to be positive 
so a < d — 1 . 

If we define 

h{t) d = f \\Su\\ + \\Sq\\ i _„_ 1 + II^H 4 _„ + ||*u|| 

#2,1 a 2,l #2,1 Lj t tS 1,\ 

+ 11*911-^ + INI £ ^-« + 2,^ (4.64) 

then collecting the previous estimates on the external terms in (4.60), there exists C v > 
(bounded when 77 G [0, 1]) such that we have for all i > (and in any case d = 2 or d > 3), 



/^) < c, e ( /V)(lky H 

v JO #2,1 



+i + IMI . d +2 + IK11 .a+i 

#2.1 #2 1 



+ ll«ll .4+1 + Ikell J+i,d)dT + KCpe 2 ^ ^\\q\\ d +2p _ a 

B 21 B e L t B 2,l 

+ (IML^ + M^i )-(\M- Ll + UHI Z1 ^ +1 -J> (4.65) 

t 2,1 t 2,1 t £ t 2,1 
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Thanks to (4.57) and (4.58), this estimate turns into: 

h(t) < CJ f h(T)(\\q\\ d +1 + \\q\\ d +2 + \\u e \\ d +1 + \\u\\ d +1 + \\q £ \\ d +1 ,d)dr 

J B 2,l B 2,l B 2,l B 2,l B £ 

+ K C^-^\\q\\ 71 .| +2 ,_ a + 4r,h(tj). (4.66) 



The term ||g||_ . d +2(3 _ Q can be estimated by interpolation (thanks again to (4.58)) if 

L t B 2,l 

1 < 2/3 — q < 2 and in this case it is less than 2C.rj. Let us collect all the conditions on 
a (see remarks 15 to 19) and f3: 

j a e]0, 1], a < |, a<d-l, 
\/3€]l,2[, ± + §</3<l + §. 

This obviously implies that d > 2, and if a is fixed in ]0, 1[ (when d = 2) or in ]0, 1] 
(when d > 3), and if we choose /3 = 1 + ^ all the conditions are satisfied. Moreover if 
T] > is such that 4r].C v < \ then the last term in 4.66 can be absorbed by the left-hand 
side, and tanks to the Gronwall lemma, for all t > 0, we have 



So (Nl .4 + i+NI d +2 +\M\ d +1 +IHI d +1 +\\le\\ d +1 d)dr 

h(t) < C v , aK e a e 



And finally, thanks to (4.58) and (4.57), we obtain that for all teR, 

L t B 2,l L t B 2,l L t B 2,l L t B 2,l 

+ \\Sq\L d a+1 d_ a + \\5q\\ d_ a+2 d_ a < C v , aK e a . (4.67) 
which ends the proof of the theorem. ■ 

Remark 20 // we had not estimated this way the terms 5G2 and 5Gs, we would write 
that 5G2 + SG3 = V(L(q £ ) — L(q)) and then due to the regularity of 8q, we could only 
obtain that \\5G2 + 5G2W . d_ a _ x < tC ri \\5q\\ . d_ a , which prevents any global in time 

L t B 2,l L T B 2,1 

estimate. 



5 Appendix 

In this section we state and prove some estimates involving the hybrid Besov spaces 
introduced in the beginning of the paper. These estimates are very close to the one 
proved in [10] (see appendix). 

Proposition 8 Let s G R, a > 0. For all q G Ba'°° D B^ 1 , we have 

h\\%s 2 i < Ikllsj.oo Ikllsji 
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Proof : This proof is classical (see [10] appendix). Prom the definitions given in the first 
section we have: 



lei 

= E (^ max ( a ,2-')||A i9 || 11 )i(^^||A„|| L2 )J 

lei v ' ' 

i i 

/ \ 2 / r.;., \ 2 



< f^2-max(a,2-)||A^|| L2 ) [^ -^—^ ^^ 
Kiel 



max^a, 2 ') 

< HgH^s.oo (5.68) 

Remark 21 More generally, if s G M, a > 0, r, r' G [l,oo] with £ + p. For all q G 

5a r n 5 a ' r , we have 

<? rjs < Wq\\ 5s, t- (/ g s r ' 
"-fc>2,l B o -Da 

Proposition 9 Let s < (i/2, t < d/2 — 1 such that s + t > 0. There exists C > such 
that for all (u,v) G 5a°° x 5a' 1 , 

ll^ll .i ( d < Cllull Ss.oo llvll St,l . 

Proof : Using the Bony decomposition, we have uv = T u v + T v u + R(u, v) (we refer to 
the introduction for the definitions of paraproduct and remainder) and we will separately 

■ s+t— 

estimate these terms in B 2 i 2 • Before that let us state the following lemma (we refer 
to [10] Proposition 5.3.) 

Lemma 4 Let a > 0, a, b G M.. Then we have 

max(a,2" a ) f 1 if a > b, 
max(a, 2~ b ) ~ \2 b ~ a if a < b. 

1- Let us begin with T u v = J2 q ez S q -\u.A q v . As for all q G Z, S q -iu.A q v has its 
frequencies localized in 2 9 C (where C is a ring) we will bound the norm by estimating 
2 q( - s+t ~ ^ \\S q -iu.A q v\\ L 2 . Using the Bernstein lemma, 

2 q ( s+t -^\\S g -iu.A q v\\ L 2 < 2«( s +*-#) I £ ||A^|| L oc J .||A^|| L2 

<2^+*-f)[ £ 2 9 'i||A g /u|| L 2 J .||A g f || L 2 (5.69) 
V?'<s-2 / 

As (u, i>) G Ba°° x 5a' 1 , there exist two nonnegative sequences c, d G / X (Z) such that 
||c||/i < 1, ||c'||;i < 1 and for all j G Z, 

2 _J,s c 

||A 1 'n||r2 < : — ^-—llnll 53,00 and ||A 1 -v||r2 < 2 _J 'c'- max(a, 2 _ - ? )|b|| s s ,i, (5.70) 

J ~ max(a,2-J) B « J _ J s « 
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so we can write that 

.d-, / max(a, 2~ q ) 



^— ' y max(a,2 9 ) " Q ^ 

Thanks to the previous lemma, as q' < q — 2 we obtain 



2 9(«+*-i)||^_ lU .A,t;|| L2 < Yl 2iq ~ q ' )iS ~^ C q ') C, q \\ U \\B^\\ V \\B^- 

\q'<q-2 

If we denote 

f 2<?(*-f) if ^ > 2 
"""jo if g < 2, 

then I ^ q ~ q ' ){s ~^c q > = (a* c)„ and as c, c' G ^(Z), a G Z°°(Z) if and only if 

\g'<<z-2 / 
s < | then (a * c).d G Z*(Z) when s < § and 

H^muII „, t d < CII U\ \ D S,00 D s il • 

II u, || . s+t _^ _ || || M ^ 

2- Let us now turn to T v u. The same arguments as previously give that for all q G Z, 

2 ^-i)||Vi-Au|| i2 < £ 2fa-^^ c ; ^ ( r a ' 2 "2 iHig.,,||t,||g., 1 . 

y max a, 2 <? 

This time, as q' < q — 2 thanks to the previous lemma, ^^("'2-9) — 2 9-13 ' and 

2^-f)||Vi-AHlL 2 < ( £ 2M(*-^ 1 ) c? ,]4i|u||^||,||^, 

\?'<9-2 / 

Similarly, if we denote 

b _r 2 ^f +i ) if q > 2 

9 ~{0 if g < 2, 

Then if t < \ - 1 we get that ||(fe * c).^^^) < ||6||joo(z) ||c||ji(z) ll^lb^z) and 

II^V^II .J-t d < Cllnll Ss.oo llf II • 
II II aS+i-j — II ll±f a II H_B a 

3- Finally, let us look at the remainder R(u, v) = A q u(A q -iv + A q v + A q +\v). Here when 
g £ Z the frequencies are only in the ball 2 q B and we estimate, for all j, 2^ s+t -^\\A j R( u, v ) 
If we perform the same computations as before (using \\uv ||^2 < IMlLira/tj/IMIi, 2 an d the 
Bernstein lemma to estimate the L°°-norm by the L 2 -norm as the freuencies are local- 
ized) , we will end with the condition s + t > | . We can get a better condition if we first 
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estimate the L 2 -norm by the Z^-norm (using the spectral localization) and then use the 
Hlder estimate. Thanks to the Bernstein lemma: 

2^ s+t -i)\\A j R(u,v)\\ L 2 < 2^ s+t -^2^\\A j R(u,v)\\ L i 

<2J(«+*) £ ||A,u|| La (||A,_i V || L2 + ||A,«|| L2 + ||A, + it;|| La ). (5.71) 

q>j-N 

As (u,v) G Ba'°° x Ba' 1 , we can use (5.70) and get: 



W\\ 5s,l 



2^+*-f)|| A , J R(^ v )|| i2 <2^ S+ *) V — £_4_||„||- ii00 x 2-9*| 
11 3 1 ' ~ ^ max a, 2-9 1 " s <* 1 

x (^2*c / 9 _ 1 max(a,2- (9 - 1) + c' q max(a, 2~ q ) + 2 - V fl+1 max(a, 2~ (c,+1) ) . (5.72) 
Thanks again to lemma 4 and (5.70), we obtain: 

2^ s+ *-f)||A^(«, V )|| i2 

<(l + 2* + 2- t )||n||^o ||t;|| 5 M £ 2tf-*>< a+t ) Cg (c',_ 1 + c' ? + c' q+1 ) 

q>j~N 

<(l + 2 t + 2- t )\\u\\^ a ,oo\\v\\^ a ,i ^ q){s+t) c q - (5.73) 

q>j-N 

Once again, if we denote 



d q = 



2 i( s +*) if j < at 
else, 



As we want (c£ * c) G Z X (Z) we now need that d G ^(Z) (d G /°°(Z) is not enough as we 
only could estimate c^_ 1 + c' 9 + d q+1 < 1 in L°° instead of L 1 ) which is true if and only 
if s + t > 0. So when s + t > we finally get: 

which ends the proof of the proposition. ■ 
Remark 22 For all q £ B^ 1 n B s 2 l = B^°° we have 

< + \\<i\\b° a 

and thanks to remark 9, when e > is small enough, for all q G Bp~ 1,s , we have 

so we can use the hybrid norms introduced in (1.13) and we will in fact use the following 
results: 
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Proposition 10 Let s£l. There exists a constant C > such that for all e > 0, and 

all q G B^ 1 n BJ,i n S| +1,s , we have 



Proposition 11 Let s < d/2, t < d/2 — 1 such that s + 1 > 0. There exists C > such 
that for ali e > 0, and all (u,v) G (B^ 1 n B^) x £* +M , 

HHL.+«-$ < C(||u||^.-l + ||«||b. )||u||Bt+l,t. 
B 2,l 
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